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ABSTRACT 

Relativistic black-hole jet sources are leading candidates for high energy (3> 
TeV) neutrino production. The relations defining (a) efficient photopion losses of 
cosmic-ray protons on target photons and (b) 77 opacity of 7 rays through that 
same target photon field imply clear multiwavelength predictions for when and 
at what energies blazars and GRBs should be most neutrino bright and 7-ray 
dim. The use of multiwavelength observations to test the standard relativistic 
jet model for these source is illustrated. 

Subject headings: gamma-rays: bursts — neutrinos — theory 



1. INTRODUCTION 

Multiwavelength observations provide important information about radiation processes 
and properties of astrophysical sources that cannot be obtained by observ ations in a narrow 
wave band. The opening of the high-energy neutrino window (for review, see lLearned fc Mannheim 
2 OOOl ) will provide a new channel of information that, in conj unction with photo n observa- 
tions, will test models of these sources. The km-scale IceCube ( lAhrens et al.ll2004l ) reaches a 
total exposure of ~ 1 km 3 -yr in 2009 and its design sensitivity in 2011. Plans are also being 
made for a northern hemisphere KM3NeT neutrino telescope in the Mediterranean Sea[j] 



Because of their rapid variability, large luminosity, and relativistic outflows, gamma- 
ray bursts (GRBs) and blazars h ave been considered as the most likely sources of ultra-high 
energy cosmic rays and neutrinos (IWaxman fc Bahcaillll997l ; IVietrilll998l ; iRachen fc Meszaros 
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19981 ; lAlvarez-Muriiz et al.ll2000f). The expected low neutrino- induced muon even t rate, even 



for the brightest 7-ray blazars (lAtovan &: Dermerll200ll ; iNeronov fc Semikozll2002l ) and GRBs 
(IDermer &: Atoyanl 120031 ; iGuetta et al.l 12004 ) , and t he increasing importance of the cosmic - 



ray induced neutrino background at lower energies (IKarle et al.l 120031 ; iGaisser et al.l Il995l ) 



means that event detection is greatly improved by choosing appropriate time windows during 
periods of highest neutrino luminosity. The important time windows for neutrino detection 
from blazars might be thought to be during 7-ray flaring states and, for GRBs, around 
the t 90 time or during some hours around the burst trigger, as suggested by the extended 
100 MeV - GeV emission for GRB 940217 (IHurley et al.lll994T) and delayed anomalous 7-ray 



emiss ion components in GRB 941017 and a few other GRBs (iGonzalez et al.ll2003l ; IGonzalez 



20051 ). because that is when energetic particle acceleration is most vigorous. But photopion 



production is enhanced in conditions of high internal photon target density, so that times of 
most favorable neutrino detection could also be argued to take place during periods of low 
7-ray flux as a result of attenuation by the dense internal photon gas. 

Here we explore this issue, considering how to use multiwavelength observations (at 
optical, X-ray, and 7-ray energies for blazars, and at X-ray and 7-ray energies for GRBs) 
to define the most favorable conditions for efficient neutrino production. Violation of these 
predictions will call into questions the underlying assumptions currently used in models of 
GRBs and blazars. 



2. Analysis 

Blazars and GRBs are widely thought to be black-hole jet sources powered ultimately by 
accretion onto a black hole, or by the spin energy of the black hole. In both cases, observations 
show that collimated outflows of highly relativistic plasma a re ejected by pr o cesses takin g 



place from compact regions. In the internal shock model (IMeszarosl 120061 ; iPiranl 120051 ) . 
collisions between faster and slower shells dissipate directed kinetic energy in the form of 
field energy and accelerated particles that radiate. After the collision, the energized shocked 
fluid shell expands on the light-crossing time scale r' b /c or longer, where r' b is the characteristic 
size of the radiating fluid element in the comoving frame, assumed spherical and isotropic 
in the comoving frame. The causality constraint implies that the size scale of the emitting 
region r' b < c8^,t var /(l + z), where the measured variability timescale t var = 10 T t T s, and 5v 
is the Doppler factor. 

Within this geometry, the relation between the measured vF v flux f e and the target pho- 
ton emissivity f{e',Q!) = dS' /dV'dt'dVl'de' is f e 5^Vyf(e',n')/d 2 L , where e = hv/m e c 2 , 
primes refer to comoving quantities, and unprimed quantities to measured values. The lumi- 
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nosity distance d>L = 10 e de cm, can be calculated for the standard ACDM universe (h = 0.72, 
Q m = 0.27, Qa = 0.73) for a source at redshift z, e'j'(e', Q') = cu' e , / 4iir' b for radiation emitted 
isotropically in the comoving frame, and u' e , = m e c 2 e' 2 n'(e') is the spectral energy density of 
the radiation field. 

We write the vF v flux as f e = f epk S(x), where S(x) is a spectral function of the variable 
x = e'/e' pk = e/e p k. Here e p k = 10- 7 is the measured photon energy (in units of m e c 2 ) of 
the peak of the uF v spectrum with peak flux f epk = 10^ ergs cm" 2 s _1 . Thus e'n'(e') = 
3^|/e pfc 'S'( a; )/( cr b 2 ^D m eC 2 e'). The rate at which protons lose energy through photohadronic 
processes is t'^J = e'n'(e ')ac, where a = 70 nb is the product of the jp photohadronic cross 
section and inelasticity (lAtoyan fe Dermerll2003l ). and the threshold condition 2^' > e' thr = 



400 relates the proper frame proton Lorentz factor 7^ and the internal photon energy. 



The target comoving photon spectral energy distribution (SED) from quasi-isotropic 
emissions, whether it is the synchrotron and synchrotron self-Compton (SSC) fields or a 
cascade radiation field, is approximated as a broken power law S(x) = x a H(l—x)+x b H(x—l) 
with vF v indices a and b (see Fig. 1; more general spectral forms can easily be treated). Here 
the Heaviside function H(x) = 1 for x > and H(x) = for x < 0. The photopion energy- 
loss rate of ultrarelativistic protons with Lorentz factor 7' interacting with photons with 
energy e' pk near the peak of the vF v spectrum is, from the proceeding considerations, 

3adlf epk (l + z) 
P * rn e c^ D tl ar e pk ' UJ 

For the model target photon spectrum with < a < 3, b < 0, 

2y b - 1 /[(l-b)(3-b)\, y»l 
C'W) " P*« { % a ~V[(l " a) (3 - a)], y < 1, < a < 1 (2) 
(a-b)/[(a _1)(1_6)], y<l, l<a<3 



(jPermerl 120071 ). where y = e' thr /2-y p e' pk = S^e' thr /2j p (l + z)e p k, and the Lorentz factor 7 P = 
E p /m p c 2 of an escaping proton as measured by a local observer is 7 P = 8d1 p - The condition 
y = 1 for protons with energy E p = E^ w interacting with photons with energy e pk implies 
that 

^ rn p c 2 5 2 D e' thr ^ 1.9 x 10 14 ^ 
p 2(1 + z)e pk (1 + ^)e pfc (keV) ' 1 } 

The radiating fluid element will expand explosively following its rapid energization 
through shell collisions, or through external shocks formed when the outflow sweeps through 
the surrounding medium. Photopion processes can be certain to be efficient — assuming 
of course that ultrarelativistic protons are accelerated in black-hole jets — if the photo- 
pion energy-loss rate p<p n , eq. (0Q), is greater than the inverse of the light travel timescale, 
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(1 + z)/8i)t V ar- An energetic cosmic ray will therefore lose a large fraction of its energy into 
electromagnetic and neutrino radiations through photopion production when the jet Doppler 
factor 



<5d < Os-k 



3&dU 



pk 



' 1 - io~ 10 - 64+ ^ 2 ^ +7? ~ r "^ //4 = 7 3d!/ 2 ( ^~ 10 ) 1 / 4 

to^pk 



(4) 



The same radiation field that functions as a target for photomeson processes is a 
source of 77 opacity. The photoabsorption optical depth for a 7-ray photon with en- 
ergy E\ in a quasi-isotropic radiation field with spectral photon density n'(e') is t. 



r' b de' <7 Tr (e / , ei)n'(e'). Using a 5-function approximation cr 77 (e / , e^) = a^e'Sie 1 
for the 77 pair production cross section <r 77 (IZdziarski &: Lightmanlll985l ) gives 



77 



pk ■ 



pk ■ 



2/^/3 

(5) 



(IDermerl 120051 ) (see also iLithwick fc Saril 1200 ll ; iBarwick et al.l 120061 ) , where 



T; 



pk 



and 



77 



pfc 



Aui e c^t v 5p 6p/j 



2<5£ 



(6) 



(7) 



1 + 

Fig. 1 compares the 5-function approximation giy e n by the a bove equations wi th accurate 
calculations using the results of iGould fc Schrederl (119671 ) (see iBrown et al.lll973L for correc- 
tions). 

At the Doppler factor 5d = 5<j, n that allows for efficient photopion production, the 77 



optical depth at photon energy e± k is, after substituting eq. 



into eq. (jSJ) 



77 



12a 



(8) 



Whenever photopion production is important, 7 rays with energies given by eq. j7j) have to 
be highly extincted by 77 processes when interacting with peak target photons with energy 
~ e p k, making it impossible to detect 7 rays at these energies. This energy is 
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The energy of protons that interact most strongly with peak target photons through pho- 
topion processes under conditions when photopion processes must be important is, from eq. 



3. Results 



Table 1 lists the important quantities derived in this paper. The quantity 5^ n is the 
jet Doppler factor where photopion losses are guaranteed to be important for protons of 
escaping energy E^ 71 . Protons with this energy undergo photopion interactions primarily 
with peak target photons with energy ~ e p k- EJ 1 is the energy of 7 rays that are attenuated 
through 77 pair production primarily by peak target photons. 

Consider the target photon variability time for the following source classes: flat spectrum 
radio quasars (FSRQs) and GRBs, both known sources of GeV radiation, and X-ray selected 
BL Lac objects (XBLs), of which over a dozen are known TeV sources. We define the 
variability time scale t var as the measured time over which the absolute flux varies by a factor 
of 2; if a source varies by N% over time At, then t var = lOOAt/iV, keeping in mind that this is 
a conservative assumption for temporal variability given that quiescent or unrelated emissions 
can add a separate slowly or nonvarying background. R-band optical and RXTE PCA (« 2 
- 60 keV) observations of 3C 279 and PKS 528+134 show that day timescale optical an d 
X-ray variability can be expected for FSRQs (IHartman et al.ll200ll ; iMukherjee et al.lll999l ). 
Day timescale optical/UV and X-ray va riability is also observed for TeV BL Lac objects 
Jpian et al.lll997i : iBlazeiowski et al.ll2005h . 



For canonical FSRQ values taken from observations of 3C 279 or PKS 0528+134, Table 
1 shows that photopion production is already important at Doppler factors of ~ 9 - 16 during 
times of day-scale optical flaring, and these optical photons effectively extinguish all 7 rays 
with energies > E^/800 1 ^ 1 -^ (Fig. 1), which would certainly include all > 100 GeV - TeV 
photons. Unfortunately, TeV telescopes have not so far been successful in detecting FSRQs, 
but monitoring of an FSRQ during an optical flare with low-energy threshold air Cherenkov 
telescope such as MAGIC would identify periods of likely neutrino emission. Photohadronic 
neutrino secondaries have energies E v « E t /20 and so would b e produced at ~ 10 17 - 10 18 
eV, providing possible sources for ANITA (IBarwick et al.l 120061 ). though outside IceCube's 
optimal energy range. 



For guaranteed importance of photohadronic production implied by day-scale X-ray 
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variability, the Doppler factors of FSRQs and TeV BL Lac objects like Mrk 421 have to be 
unexpectedly small, < 3. If the X-ray flaring timescale of FSRQs were hourly rather than 
daily, then 5^ would more nearly correspond to Doppler factors ~ 5 - 10 inferred from uni- 



fication studies and sup erluminal motion observations of blazars (jVermeulen &: Cohenlll994 



Urry fc Padovanilll995l ). During episodes of highly variable X-ray flux, such sources should be 
invisible to GLAST, and ^> TeV neutrinos should be created if FSRQs are sources of ultra- 
high energy cosmic rays. For the XBL estimate, a 15 minute X-ray flaring timescale has 
already been assumed. Thus FSRQs are more likely than BL Lac objects to be high-energy 
neutrino sources for IceCube, which also fo llows if, as is likely, the external radiation field 
plays a stro ng role in neutrino production feednarek Sz Protheroe 1999 ; Atoyan Sz Dermer 
200lLl2003h . 



The outcome of this analysis to identify periods of high-energy neutrino production is 
best for bright GRBs with peak fluxes of ~ 1CT 6 ergs cm -2 s™ 1 and peak photon energy in 
the range 50 keV - 0.5 MeV that show < 1 s spikes of emission. The b ulk factors, ~ 100, are 



consistent with widely considered outflow speeds in GRBs (see, e.g.. iRazzaque et al.l 12004 



who also treat 77 attenuation in GRBs). Perhaps 100 MeV photons could be observed, but 
the GLAST LAT should see no > GeV photons if <5d ^ $4™, which is the most favorable time 
for detecting 100 TeV - PeV neutrinos and is at an optimal energy for detection with km- 
scale neutrino telescopes. Bright X-ray flares with durations ~ 10 2 s observed hundreds to 



thousands of seconds after the GRB trigger, like tho se discovered with Swift (jO'Brien et al. 



20061 ; iBurrows et al.ll2007l ; iMurase &; Nagatakill2006l ). with blast wave Doppler factors ~ 50, 
are also promising times to look for neutrinos and a 7-ray spectrum attenuated above ~ 100 
GeV 7 rays. 



The condition r 77 (ei) < 1 for detected 7 rays with energy e\ gives a minimum Doppler 



factor 



°D 



4m e cH v e£ k 



l/(6-2A) 



(11) 



(fPermerl 120051 ). with A = b if e x < ef , and A = a if e x > ef . If 5% in > 5^, we should 
not expect GRBs to be neutrino bright. Synchro-Compton analysis of high-quality radio 
and gamma-ray blazar SED data with resolved VLBI cores and self-absorption frequencies 
give the Doppler factor directly, provided the synchrotron self-Compton component can be 
identified and separated from external Compton and photohadronic emission components. 
Should high-energy neutrinos be detectected when the Doppler factor inferred from these 
tests is greater than 6^, then this would call into question our understanding of the structure 
of black-hole jets, for example, the assumption of isotropy of target photon distributions in 
the comoving jet frame. 
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4. Summary and Conclusions 

We have presented a detailed treatment of combined photomeson and gamma-ray opac- 
ity, which is a crucial problem that unites the neutrino particle physic s and the electro- 



magn etic worlds. This problem has been numerically treated for GRBs (IDermer fc Atoyan 



20031 ). but no detailed analytical treatment is present in the scientific literature. GLAST 
observations will reveal if 7-ray spectra of FSRQ blazars and GRBs show evidence for strong 
7-ray absorption during periods of variable target photon emissions, signifying favorable 
conditions for high-energy neutrino production. 

To illustrate the use of these results, suppose that a blazar or GRB with measured 
redshift z is monitored at optical, X-ray, or soft 7-ray energies, giving a light curve f e (t) at 
photon energy e(t). The structure of the light curve implies t var (t). From these observables, 
the Doppler factor 5</> n (t) for guaranteed importance of photomeson production is derived 
from eq. PJ. If high-energy neutrinos are detected, then the source must be opaque at 7-ray 
energies given by eq. (jUJ). If 7 rays are detected at these energies, then the basic relativistic 
jet model must be wrong. Suppose instead that 7 rays at some energy e x are detected. In 
this case, the minimum Doppler factor 5™ n (t)can be inferred from eq. (11 ip to define times 
when these sources can and cannot be neutrino bright. 

Times and locations of bright, variable MeV 7-ray and extincted GeV fluxes from 
GRBs can be done exclusively with the GLAST GBM and LAT, whereas other tests for 
7-ray/multiwavelength correlations giving the most favorable times for high-energy neu- 
trino detection in blazars require collaboration and coordination of separate facilities. The 
necessary organization is already underway between GLAST and the ground-based 7-ray 
telescopes, e.g., HESS and VERITAS, but blazar observations with, e.g., Swift, Suzaku, and 
RXTE correlated with GLAST, AGILE, and ground-based high-energy 7-ray telescopes will 
be crucial for neutrino discovery science and testing models of relativistic jet sources. 
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Fig. 1. — Target photon SED, 77 opacity r 77 (ei), and normalized photopion energy-loss rate 
t'dm / P<t>-* are shown as a function of photon or neutrino energy for parameters of (a) optically 
flaring FSRQ and (b) rapidly varying prompt GRB emission, using parameters from Table 
1. The target photon SED is approximated as a broken power law with uF v flux peaking 
at e p k- When photopion processes are certain to be important for protons with energy 
that interact with peak target photons with energy fa e p k, then the 77 opacity of 7 rays with 
energy E^ 7 is ~ 800. The 77 opacity is less than unity at photon energies 
Heav y and light curves for r 77 (ei) ar e accurate nu merical integrations (IGould fe Schreder 



19671 ) and 5-function approximations (jDermerll2005l ). respectively. 



Table 1. Doppler Factor 5^, 7-Ray Photon Energy E^ 1 , and Cosmic Ray Energy E^ n 
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30 x 10 3 


1.6 x 10 19 
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5 x 10 16 


X— ray flare 




-9 


2 
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290 


1.6 x 10 17 



Sources at z = 2 except for XBLs, at z ~ 0.08, di = 10 27 cm. 



